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The new approach to the microscopic description of the phase transitions starting from the only 
first principles was developed on an example of the transition normal metal-superconductor. This 
means mathematically, that the free energy is calculated in the range of temperatures, which in- 
cludes a point of pase transition, without introducing any artificial parameters similar to an order 
parameter, but only starting from microscopic parameters of Hamiltonian. Moreover the theorems 
about connection of a vacuum amplitude with thermodynamics potentials are realized. The func- 
tional of a superconductor's free energy in a magnetic field was obtained with help the developed 
method. The obtained functional is generalization of Ginzburg-Landau functional for the case of 
arbitrary value of a gap, arbitrary spatial inhomogeneities and nonlocal magnetic response. The 
explicit expressions for the extremals of this functional were obtained in the low-temperature limit 
and the high-temperature limit at the condition of slowness of gap's changes. 
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I. FORMULATION OF THE PROBLEM 



Exact microscopic description of a phase transition is opened problem of modern theoretical physics. The problem 
can be solved for several simplest models only, but different phenomcnological approaches exist for the rest cases. The 
essence of the matter lies in the following. Basic problem of statistical mechanics is calculation of a partition function 
(of grand canonical ensemble in a total case) Z or calculation of a density matrix p: 

Z = Sp(e-^-^)=e-P n , p=L e -^ S -^\ (1) 

where H = Hq + V is full Hamiltonian of a system, N is the particle operator, p is the chemical potential, f2 is the 
grand thermodynamics potential. Replacement of Hamiltonian H in canonical ensemble by Hamiltonian H — pN in 
grand canonical ensemble leads to the shift of reference of particle's energy from zero to Fermi surface: e{kp) = 0. 
Hence, the potential plays a part of Hclmholtz free energy with a reference of particle's energy from Fermi surface. 
Therefore we shall call the grand thermodynamics potential by free energy for brevity. Let's transform the partition 
function fTJ) to a form: 

Z = Sp (e-W'-^ePV'-^e-KB-i^) = Sp (z p U(p)) = Z R(/3), (2) 
where Zo is the partition function for a system of nonintcracting particles, 

U { P )=e +m -,N )e -m-,N) =gizil: I" dTl ... f dr n f{H I (r 1 )...H I {r n )} (3) 

n=0 H ' J° ^° 

is the evolution operator in the interaction representation (it describes evolution of the system in imaginary time 

it — » r, T is the ordering operator in time), Hi(j3) — e Jr ^ tyHo ~^ N ^Ve~^ tyHa ^^ N ^ > is the interaction operator of particles 
in interaction representation, 

R(fl) = (U(fi)) = Sp (poUtfj) = —T J Q dT i ■■■J Q dT nS P (pbf • ••£/(«}) (4) 
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is the vacuum amplitude of the system. The averaging ()o = Sp(fio . . .) is done over ensemble of noninteracting 
particles. It is necessary to note, that the vacuum amplitude describes dynamics of the system in the multiparticle 
state $o under the influence of the internal interaction, unlike Green function describing dynamics of a particle in 
one-particle state 4>k under the influence of its interaction with other particles. 

The partition function Zq can be found exactly for any system. If particles interact, then the situation becomes 
complicated essentially. The only several models can be solved exactly @, E|, but in the rest cases the values (p} 
can not be found exactly and it is necessary to develop a perturbation theory. Solution of the basic problem of the 
statistical mechanics reduces to the calculation of a transition amplitude " vacuum- vacuum" Q. In order to formulate 
our problem let's consider the case of zero temperature T = 0. In this case the time is real, and the vacuum amplitude 
is determined by the following: 

R(t) = (<S> \U(t-t )\<5> ) to=0 = (<f> \U(t)\<f> Q y Wot 

= J2 J dtl ---J dtn • • ■ l $ o) , (5) 

where Hi(t) = e + lt ( H o- ^ N )y e -it(H - ^n) - m Q p era tor of particles' interaction in the interaction representation, Wo is 
ground state energy of a system without interaction. The averaging ()o = (<I>o | ■ ■ • I'&o) is done over the ensemble of 
noninteracting particles. Knowing a vacuum amplitude we can calculate the ground state energy of the system of the 
interacting particles Eq using the theorem [l|: 

E = W + lim i—kiR(t), (6) 

t->oo(l-j<5) dt 

where S is infinitely small value, but oo • S — > oo. The theorem is correct if the limit transition 

is possible. This transition is possible at the condition when the symmetries of ground state of the system with 
interaction l^o) and without interaction |<i>o) are identical: 

($o|*o)^0. (8) 

The expressions ([7]) and ([5]) mean: l)potential of interaction is being switched slowly in the system in the ground 
state without interaction, 2) the ground state of the system with interaction is being obtained by continuous way 
from the ground state without interaction while the switching of the interaction (adiabatic hypothesis). For nonzero 
temperature the analog of the theorem ([6]) has a form: 

U = -±UZo-±UW), (9) 
where U is the internal energy of a system. Moreover, we can calculate the free energy: 

n = -~lnZ -~lnR(/3). (10) 

If the wave functions are orthogonal (^ol^o) = - the adiabatic hypothesis is not valid, then the symmetries of the 
system with interaction and without one is different. This means, that an initial system without interaction suffers 
phase transition stipulated by the interaction. A nonfulfilment of the adiabatic hypothesis means nonfulfilment of the 
theorem (|6]). For the system with the broken symmetry we can calculate a vacuum amplitude on the free propagators 
Go(k, t) and use the formula ©• However we shall find a wrong ground state energy E . This means that a state 
exists with more low energy than the found value. Moreover, in consequence of the breakdown of the condition © a 
system becomes unstable: T- matrix, which determines vacuum amplitude R(t), one-particle propagator G(k,i) and 
two-particle propagator if(ki,k2,t) in stair approximation, has the form: 

T(t) =ce- at +c'e +Qt . (11) 

The value of T is increasing infinitely at t — > oo 0, Q , it means an instability of the system. Similar instability was 
observed experimentally as the process of formation of a charge density wave in TbTe3 @. 
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From the aforesaid we can see, that perturbation theory making possible to calculate the vacuum amplitude (ther- 
modynamics function U , Q and so on) in ranges including a point of phase transition (for example, temperature Tq) 
doesn't exist. In other words, it is not possible to find R{f3) in the system with condensed phase starting from the first 
principles. However the phenomenological approach exists for calculation of tt. It assumes, that at the temperature 
T < Tc the condensed phase is characterized by some order parameter r\ exists. Near a transition point for type 
II phase transition or near a point of overcooling for type I phase transition the parameter 77 is small, and the free 
energy can be represented in a form of Landau expansion: 

Q(T, V, fi,h) = Q + J d 3 r(arf + brf + g(Vr]) 2 - r)h), (12) 

where h is external field. The equilibrium value of order parameter 77 is determined by an extremal of the functional 
(IT2"j) : jl\ n =rj = 0. However, the expansion (fT2j) is correct in the range (Tc — T)/Tc -C 1 only, hence rj can be found 
in this range only. 

In order to obtain rj at any temperature the concept of quasi- averages (anomalous averages) is introduced. In 
the paper Q generalization of the method of quasi-averages has been represented in terms of Green function - 
Nambu-Gor'kov formalism @, It is postulated, that in the condensed phase, in addition to normal propagators 
G(k, er, i' — t) = — i(^ \T{C ka (t')C^ (t)}\^ ), anomalous propagators F exist, for example: 



Ferromagnetic 
Solid, liquid 
Superconductor with singlet pairing 



F fer = -<(*o|T{(7 fct (f)t3+ (t)}|* ) 
F sol = -i(Vo\T{C k+(la (t')C+ a (t)}\M>c 
F sup = - t (^o\T{C_ ki (t')C kt (t)}\^ ) 



(13) 



They are proportional to the according order parameters. The anomalous propagators F can not be obtained by 
summation of diagrams consisting of normal propagators G only. This fact is result of different symmetries of a 
perturbed state and an unperturbed state. 

For construction of the self-consistent perturbation theory the "sourse term" Hs is introduced in Hamiltonian 
instead of the interaction operator V. The sourse term is induced by the order parameter (for ferromagnetic - 
magnetic field orientating spins, for superconductor - a sourse of Cooper pairs). The sourse term inducts specification 
structure and changes symmetry of a system: $0 such that (^qI^q) 7^ 0. This means that a system has an 

internal long range field H, generated by the sourse Hs, and according order parameter r\, which is function of H and 
depends on the interaction constant A and temperature T: rj = rjxr(H). In turn, H is function of 77: H = H\t(v)- 
Hence, these two equations can be combined: 



r) = ri\ T (HxT(ri)) 



(14) 



and they can be solved in 77. This means, that order parameter is determined in self consistent way. In Green function 
formalism this fact has the following form: the free matrix propagator Go (anomalous part is zero F = 0) and the 
dressed propagator G (with normal G part and anomalous F parts) obey Dyson equation 0, [E| : 



G- 1 = Go 1 - S(G) 



(15) 



and what's more the mass operator S(G) is determined by self consistent way. This means, that elements of diagrams 
for the mass operator are dressed propagators G, which contain the sought mass operator. Artificiality of the 
introduction of the sourse H$ lies in the fact that the expression for a mass operator is postulated. Moreover, X 
can not be obtained by summation of the diagram consisting of free propagators Go only, that is S(Go) = 0. For 
example, S has the form for a superconductor: 



E(w, p) = [1 - Z(w, p)]wl + x(w, p)r 3 + ip(w, p)n + <p(w, p)r 2 , 



(16) 



where Z is some coefficient, the field \ determines a shift of chemical potential /1 at the transition to the superconduc- 
tive state, the fields ip and <p play a part of order parameter in superconductivity - a gap, and specter of excitations 
is represented by them E 2 (p) = e 2 (p) + <p 2 (p) + ip 2 (p); 1, n are unit matrix and Pauli matrixes. The dependence of 
all fields on frequency u> takes into account a delay and a damping of quasi-particles. In Hartree-Fock approximation 
we have Z = 1 and S doesn't depend on uj. We can suppose tp = 0, x = and denote ip(p) = A(p), then 



(17) 
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where A < is the interaction constant, and the summation is done over Fermy frequencies uj n = (2n+ l)irT. We can 
see, that the order parameter is determined by the anomalous propagator F, and the equations (jTTJ) are equations of 
self consistency for order parameter A as specific case of the total equation (fT4| . 

The described approach is phenomenological too, because the anomalous propagators and corresponding order 
parameters are introduced to the theory from elsewhere. This means that we select the required states from all 
possible state artificially. In this sense this approach likes Landau approach (fl~2|) . The sourse of Cooper pairs for 
superconductor or the magnetic field H = (S z ) for ferromagnetic (Jo is the exchange integral, g is the gyromagnetic 
relation, \xb is Bohr magneton, (S z ) is the average projection of spin onto axis z) can not be interpreted as real field. 
So, for iron the real internal magnetic field is ~ 10 3 oersted, but for the ordering of spins it is necessary the effective 
magnetic field H ~ 10 6 oersted. Thus, self magnetization has nonmagnetic nature evidently. 

As it has been noted in Q , in spite of all progresses reached in description of phase transitions and in calculation 
of main characteristics of a system in critical region, the basic problem of phase transitions is not solved: calculation 
explicit expressions for thermodynamical functions of a system in ranges, which include a point of phase transition 
as function of temperature, external fields and microscopic parameters of Hamiltonian. It is necessary to have a 
description of phase transitions on microscopic level. This presupposes a direct calculation of free energy 0(T, V, h) 
from first principles, but we must not construct it. 

At the present moment the two approaches can be separated for solution of the formulated problem. In the papers 
[9MTT| a partition function of Ising model is calculated by the method of collective variables, which are oscillation 
modes of spin moment. Then a partition function can be written via these variables as a functional. Investigation of 
Euler-Lagrange equations shows, that among a set of collective variables the variable connected with order parameter 
exists. In the paper [l2j a partition function is calculated by the saddle point method for classical system with short 
range attraction and repulsion between particles. It has been shown, that the free energy can be represented by the 
form which is analogous to Landau expansion where the saddle point plays a role the order parameter. The saddle 
point method has the sense as method of separation of states giving largest contribution in a partition function and it 
is equivalent to the mean field method. With the help of the saddle point method the processes of cluster formation 
can be investigated. The cluster formation can be considered as a phase transition from spatially homogeneous 
distribution to spatially inhomogeneous distribution. In the papers [13l4l6l | both short range potentials and long 
range potentials were considered. Critical temperature and critical concentration, when clusters form in a system, 
dependence of size of a cluster on temperature have been obtained. 

One of the important applications of the microscopic theory of phase transitions is description of thermodynamics 
and electrodynamics of superconductors. It is necessary to know the functional of free energy Q(/3, A, A), where A 
is potential of magnetic field. Then the equations = 0, j^- = will describe equilibrium states of the condensed 
phase and normal phase. Two basic methods for obtaining of the sought equations exist. First of them is joint solution 
of Gor'kov equations and the equation of self consistency (fTTj) . At the temperature T — > Tq, when A/Tc — > 0, the 
solution of these equation can be represented in the form of series in degrees of A. Moreover, magnetic penetration 
depth is bigger then Pippard coherent length Iq, hence the potential A changes a little on a coherent length. As a 
result we have the well-known Ginzburg-Landau equation. 

Another method, proposed in [T^ |. is the direct calculation of a vacuum amplitude R(f3). The concept lies in the 
fact that we consider electrons in a normal metal propagating in random "field" of thermodynamic fluctuations of 
order parameter Aq, where q is small wave- vector. The operator of the interaction of electrons with the fluctuations 
can be written as: 



where p ± = p ± q/2. A correction to the thermodynamics potential from any interaction is represented via the 
vacuum amplitude R: 



Then using Wick theorem we can represent Q via the propagators. For the correction of second order we have: 




(18) 



P 



Aft = -T\nR{[3) 



T[R(f3) - 1]. 



(19) 




1/T 



1/T 





1/T 



1/T 



(20) 
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The correction Af2 is represented via the free propagators Go of normal state only - we consider normal metal at 
T > T c) where the fluctuation sourse of Cooper pair (TT8)) acts. As a result we have Landau expansion: 



where a(T) oc (T — T), b, 7 are expansion coefficients. 

In our opinion, this approach is not successively microscopic, because the artificial element is used - the external 
sourse of Cooper pairs (fT8|) . that implies some seed order parameter. As for calculation of the correction Af2 the free 
propagators Go of normal phase is used only, the condensed phase is considered as fluctuations against the background 
of normal phase. This means, that we can obtain the limit f2(T — > T) only. 

Ginzburg-Landau equations are correct for description of thermodynamics and electrodynamics of a superconductor 
at the following restrictions: 

1. The gap is much less than critical temperature. Then the parameter A(r, T)/Tc <C 1 can be expansion param- 
eter. This means, that the equations are correct in the range T Tc or in the range H — ¥ Hc2 (intensity of 
magnetic field is near the second critical magnetic field Hci)- 

2. A(r, T) changes slowly on the coherent length l(T), which is size of a Cooper pair. 

3. Magnetic field H(r) = rotA(r) changes slowly on the coherent length, that is the magnetic penetration depth 
is A(T) ^> 1(0). This means, that electrodynamics of a superconductor is local. 

In the papers [ill (3 the equations has been proposed, where the first restriction is absent. These equations were 
obtained from Gor'kov equations and they are the generalization of Ginzburg-Landau equations for the case of 
arbitrary value of A(r, T)/T. However spatial inhomogeneities are slow and electrodynamics is local. 

Our aim is the description of phase transitions on microscopic level starting from the first principles only. Mathe- 
matically this means to develop a method of calculation of the partition function (TTJ (the free energy f2) in ranges of 
temperatures and parameters of interaction, which include a point of pase transition, without introducing any artificial 
parameters of type of order parameter rj and sourses of ordering H$, but starting from microscopic parameters of 
Hamiltonian H = Hq + V only. The theorems (|6I9I10[) about connection of a vacuum amplitude with thermodynamics 
potentials must be realized. Thus, in microscopic theory of phase transitions the equation of self-consistency (|14l) . the 
anomalous propagators F(p, uS] (|I3|) and Landau functional (|I2[) must be deduced, but they must be not postulated. 
To solve the problem means to develop the perturbation theory for the vacuum amplitudes R(t) and R(/3), which is 
correct both normal phase and condensed phase. Phase transition normal metal - superconductor has been considered 
as an example. Since in Nambu-Gor'kov formalism (the method of anomalous propagators) any phase transition can 
be described then our method can be generalized to the rest transitions (ferromagnetism, waves of charge and 
spin density, crystallization and so on). Further we can apply the developed method in order to calculate free energy 
of a superconductor at arbitrary temperatures, spatial inhomogeneities, magnetic fields and currents, moreover with 
nonlocal magnetic response. The obtained expression will be the generalization of Ginzburg-Landau functional (|I2[) 
in above mentioned sense. 

In the section|TI]we consider the instability of normal Fermy system at the switching of attraction between particles. 
Mathematically this is expressed in the fact that a two-particle propagator K, calculated on free one-particle propaga- 
tors in stair approximation, has a pole a which doesn't belong to a free propagator Kq and situated on complex axis. 
This means a presence of bound states of particles in the system (with the binding energy \a\) and evolution of the 
system in time as ~ e'"'*. On the other hand, a two-particle propagator has a pole at the energy of the bound state if 
the bound state of two isolated particles exists. The residue in this pole is product of Bethe-Solpiter amplitudes 7777 + 
- the amplitudes of pairing. 

In the section IIIII we generalizes the result of two-particles problem to a many-particle system. It follows from the 
identity principle, that amplitudes of pairing are determined by dynamics of all particles of the system, and observed 
values of the amplitudes of pairing are the result of the averaging over a system. Thus, the collective (condensate) 
of pairs exists. In order to find an one-particle propagator Gs and to generalize the two-particle problem to a 
many-particle case we proposed the method of an uncoupling of correlations. The method considers interaction of a 
additional fermion with fluctuations of pairing (formation and decay of the pairs). As a result of such interaction, 
the law of quasi-particles' dispersion changes as e(k) — > y/ e 2 (k) + AA+, where A and A + are amplitudes of pairing 
playing a role of a gap and they are analog of Bethe-Solpiter amplitudes in two-particle problem. Gor'kov equations 
and the existence of anomalous propagators F and F + follow from Dyson equation for the described above process. 
This fact means breakdown of a global gauge symmetry, namely number of particles is not conserved in the course of a 
presence of a pairs' condensate. After calculation of particles' interaction with fluctuations of pairing all characteristics 
of a system must be calculated over the new vacuum with broken symmetry. 




q 



(21) 
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For calculation of observed values of the amplitude of pairing A and A + it is necessary to know a vacuum amplitude 
R(t) of a system. The vacuum amplitude must be calculated over the new vacuum with broken symmetry. This gives 
possibility to use the theorem about connection of a vacuum amplitude with ground state energy of a system. In the 
section ITVl the method of uncoupling of correlations is proposed. The method allows to represent a vacuum amplitude 
via anomalous propagators F and F + . Thus we obtain a functional of ground state energy over the fields A and A+. 
An extremal of the obtained functional is the equation of self consistency for the parameter A in Nambu-Gor'kov 
formalism. This means, that order parameter is averaged Bethe-Solpiter amplitude over all system. In the section IVl 
we generalize the results of two previous sections for the case of nonzero temperature. Using the method of uncoupling 
of correlations we calculate a vacuum amplitude R{0) and a functional of free energy f2(A, T) over the fields A and 
A + . In a high-temperature limit T — > Tq the expansion of free energy in powers of |A| has a form of Landau 
expansion. This fact proves, that the averaged over a system the amplitudes of pairing A and A + have properties, 
which are analogous to the properties of an order parameter in a phenomcnological theory. In the section IVII we 
consider the case of a pairing with nonzero momentum of a pair's center of mass. 

In the sections fVIII and rVIIII using the developed method of microscopic description of a phase transition we obtained 
the functional of free energy of a spatial inhomogeneous superconductor in magnetic field. The functional generalizes 
Ginzburg-Landau functional in cases of arbitrary temperatures, arbitrary spatial inhomogeneities and nonlocality of 
magnetic response. The free energy has a form of Ginzburg-Landau functional in the high-temperature limit at the 
condition of slowness of gap's change in space. Corresponding equations of superconductor's state demonstrate the 
nonlinear connection between the current and the magnetic field. In the low-temperature limit in the case of weak 
field H <C He at the condition of slowness of gap's change the nonlocal connection between the current and the field 
appears, which is the long-wave limit of Pippard law. The last fact proves the nonlocality of the obtained functional 
of free energy. 



II. THE INSTABILITY OF NORMAL STATE AND TWO-PARTICLE DYNAMICS 



Let we have a system from TV noninteracting fermions in volume V . In ideal Fermy gas propagation of a particle with 
momentum k, energy e uj?(|k| — kp) counted off Fermy surface (we are using system of units, where K = ks = 1) 
and spin a is described by the free propagator: 

' "WC^WCJ (0)|*o>, *>ol_ ^ w 

i<*o|C+ CT (0)C k >)|$o>, ^ n >-- l6tAoe +l6 - tBoe ' 



GoCM) = 

G(k,i) = J ^G(k, w )e— * ==► G (k,w) = 




uj — e{k) lj 2 — e 2 to — \e\ to ■ 



(22) 



Ao4fi + oV b = Ui-^), e t = {]: *>U, (23) 



where 

Cjj a (t) and Cj^ (t) are creation and annihilation operators in Heisenberg representation. Now let attractive force 
acts between particles. The force is described by matrix element of interaction potential: 

(l,-l\V\k,-k)=Xw lWk <0, = ^)>^}» ( 24 ) 

moreover interacting particles have opposite spins. This model potential is correct if a range of interaction tq is much 
smaller than average distance between particles: Vq ^/ N/V ~ r^kp -C 1, that means "slowness" of collisions. In turn, 
account of Fermy statistic, in the limit of "slow" collisions particles can scatter with opposite spins only. 

The mass operator E(Go) in stair approximation is determined by so called T-matrix (amplitude of scattering), 
which is solution of the equation represented by the diagram in FigQ] For the case q = T-matrix has a form: 

r _ \w k ,w k _ \w k ,w k , 

d 3 k dw„ ^ '■• n ~ f n ~ . . -.. A,.a ■ y Zb > 



1 



iXf ^^w k G (k,a-io),Go(k,cj) l-|A|^f In 



This expression has a pole in the point ag. In its neighborhood the expression has a view: 

r(Q -> od) = -^-^L al(X 0) = -4u£exp {-^j } • (^(i) 
mkp a. ±t|ao| L mkp\\\ 
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q-k q-k' q-k q-k" q-k" q-k 

G o( k ' w ) a-co a-ffl' a-co a-co a-co' G (q-k",a-co) a -co 




Figure 1: The mass operator £ expressed via F(q, a)-matrix. Bethe-Solpiter equation for F(q, «)-matrix. 



If with the help of Fourier transformation to pass from ^-representation to t-representation, then we shall have the 
expression (|11[) . which increases infinitely at t — > oo. Hence, according to the first diagram in Fig[TJ we have the same 
bad behavior of the mass operator E. The presence of a pole in imaginary axis means instability of a system. In our 
opinion, this result can be interpreted as follows: in a system because of the interaction (|24l) the strong fluctuations 
exists, but the free propagator Go doesn't consider these perturbations. In turn, it leads to the instable solution of 
Dyson equation. This means, that besides the interaction between particles V(k) the interaction of particles with 
aforesaid fluctuations exists too. For stability of the solution, dressed propagators Gs, which considers scattering on 
the fluctuations, must be in the equation for the mass operator E and scattering amplitude T in Fig. (|T[) instead the 
free propagators Go- 

Let's determine a two-particle propagator by the expression: 

K(x 1 ,x 2 ;x 3 ,x 4 ) = (^ \TC(x l )C(x 2 )C + (x 3 )C + (x 4 )\^ ) 

K (x 1 ,x 2 ;x 3 ,x 4 ) = (^ \Td(x 1 )C(x 2 )C + (x 3 )C + {x 4 )\^ ) = G (x l ,x 3 )G (x 2 ,x i ), (27) 

where K is a propagator in a system with interaction, K is a free two-particle propagator and it is equal to a 
product of one-particle propagators Go, x = (£, t). The free and dressed two-particle propagators are connected by 
Bethe-Solpiter equation (in an operator view): 

K = K + K Q ivK = K + K iTK => V = v + ivK T, (28) 

where v{x\, x 2 ; X 3 , x\) = V(r-L—r 2 )5(xi—x 2 )5(x 2 — X4)5(ti — t 2 )- The last equation in for the scattering amplitude 
T is represented graphically in FigJlJ We can see from the equation ([28]) . that a p resence of the pole in T means a 
presence the same pole in the two-particle propagator K. As it is well known |20L [22| - [25| . a two-particle propagator 
has a pole structure at the values of energy corresponding to a bound state. The pole |ao| (|26[) doesn't belong to the 
free propagator Kq, but it appears as a result of the attraction ([24]) A < 0. Hence, the pole means a presence of bound 
states of two particles in a system with the binding energy E s sa |a |- 

For investigation of the bound state and calculation of of particles' interaction with above described fluctuations 
let's consider the problem of two particles at first. Previously we considered dynamics of two selected particles being 
in a field of the rest particles of a system. Therefore the propagator K determined dynamics of all particles of a 
system in such approximation. Now we shall consider the system consisting from two particles only being in the state 
Q s with the energy E s . Interaction between them is described by the potential V(k) = V\ 2 : 

{H 1 +H 2 + V 12 )® s =E S $ S . (29) 
Let's determine a propagator for two particles as kernel of the integral operator finding $(£) known <I>i(i'): 

*(£i,6,t > * / ) = - / K(^^ 2 ,t^[,e 2 ,t')^[,e 2 ,t')d^de 2 . m 

Then the two-particle propagator can be written in Fourier representation as (20| : 



(31) 
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If the bound state is among states s, then the pole of the function K(^±,^2,E) corresponds to the bound state s at a 
real E, where E is equal to energy of the bound state E s . 

The residue in a pole E = E s is $ s (£i, £2)$s(£i, £2)- As ^0 hasn't a pole in E = E s , where E B is the energy of 
the bound state, then a pole of K means a presence of a pole at T according to the equation (|28p. For the function 
K(xi 1 xi\ x'±, X2) 1 with different times the expression can be written: 

^ (€1, £2, n, £1,42^1^) = 2^ E-E + i-y ' ' 

s s 

where we denoted that 



ti — t<i = r l5 t\ — t 2 = t[, t\ + t-i = 2t, t[ + t 2 — 2t' . 

This expression is the formula (pTLj) at Ti = t{ = 0. 

The residue II S (£1 , £2 , T\ ; £J , £ 2 > T \ ) f° r the bound state can be written in multiplicative form as in the case t± = 
r{=0: 

na = v»fa,&,n)rit (33) 



The values r\ and 77 + are Bethe-Solpiter amplitudes |22| - |25| . They are connected with the wave functions <1> S as (in 
momentum representation £ = k): $ s (ki,k2) = J 7? s (ki, k2, ejsp- In the equation for if (|28)l . i\To can be neglected 
near a pole corresponding to a bound state. Resulting homogeneous equation has the solution moreover the 

function t? s (£i, ^2> Ti) satisfies the equation 

/de' dk 
T/(fc)77 s (ki+q,k 2 -q,e')7^77r^- ( 34 ) 
Z7T (ZTT )° 

The term in the sum (|32[) . corresponding to the bound state -Es, can be represented by the diagram in Fig 121 where 



+ 

K ^ v » k+q 



1l(k,k 2 ) r) (k+q,k 2 -q) 



d- - a: 



k,-q 



Figure 2: The two-particle propagator for isolated pair of particles in neighborhood of a pole corresponding to the bound state 
E s . 

the dotted line means the multiplier +il - propagation of two particles in bound state, and the blocks are 

i] s (ki , k2 , e) and 77+ (ki + q, k2 q, e) - transition amplitudes in bound state and back. The diagram in Fig(3] can be 
interpreted by the following way. Two particles with momentums ki and k2 form a bound state with energy E s with 
amplitude ?7s(ki, k 2 ). Further, the two particles propagate together. Then the bound state can decay with amplitude 
r/f (ki + q, k2 — q) . As a result the two free particles appear with momentums ki + q and k2 — q. 



III. THE UNCOUPLING OF CORRELATIONS AND A MULTIPARTICLE DYNAMICS. 

In the previous section we considered dynamics of two isolated particles. Now we have to generalize the obtained 
results to the multi-particle case - propagation of two interacting particles in a system of identical fermions. This 
situation differs from the previous case by the following conditions: 

1. Each pair of fermions is in field of all the rest particles. 

2. All particles of a system are identical. Moreover, the average size of a pair l ~ l/y/\ocn\2m » yV/N is more 
big than average distant between particles, that means the wave packages of pairs overlap strongly. 
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Mathematically this means, that the amplitudes r\ and r) + are not solution of the equation (|34|) . which is correct for 
isolated pair only. Now the amplitudes are determined by dynamics of all particles of the system, and their observed 
value is result of an averaging over a system. Pairs in such system are effective, namely two fermions having formed 
a bound state with an amplitude ?y s (ki, lt2, e) (FigJ5]) are not fixed pair: one from partners in a pair can leave the 
bound state with a fermion from another pair with amplitude 77+ (ki + q, k2 — q, e). Thus, the collective of pairs 
(condensate) exists. 

Let's consider the two-particle propagator Ke^e s represented in Fig[21 As it has been noted earlier, fermions with 
opposite momentums and and opposite spins form a pair. Let's suppose, that corresponding amplitudes of pairing 
don't depend on time. In order to obtain an one-particle propagator Gs we shall use the method of uncoupling 
of correlations considered in the Appendix [A] and we shall be acting analogously to FigfTT] The procedure of an 
uncoupling is represented in Fig(3] We connect the entering line and the outgoing line corresponding to particles with 
momentum — k and energy parameter —uj each. As a result we have the intermediate propagator Go- Since partners 
in each pair are not fixed in consequence of the identity principle and strong intersection of wave packages of pairs, 
that is a condensate of pairs exists, then it is necessary to cut the dotted line - the propagator of a pair E _^ +i . Then 
the points of a joining of the dotted lines correspond to interaction with the effective field (in accordance with the 
rules of diagram technics in the Appendix EJ. The fluctuation of pairing play a role of the above mentioned effective 
field. The amplitudes of such interaction we denote as — iA(k, — k) and iA + (k, — k). These values correspond to the 
amplitudes of the two-particle problem r7 s (k, — k) and f?^(k, — k) to the extent that their observed values is result of 
averaging over a system (r] s (k., — k)) ~ A (77+ (k, — k)) ~ A+ in consequence of statistical correlations between pairs 
and they are determined by dynamics of all system's particles. 




Figure 3: The procedure of uncoupling of correlations for a two-particle propagator of a pair being in a field of all rest fermion 
of a system. The result of the uncoupling is the dressed one-particle propagator Gs(k, uj), as a consequence of interaction of a 
free fermion with fluctuations of pairing. 

The result of the procedure of uncoupling of correlations means the follows. Let an additional particle with 
momentum k, uj propagates through a system of identical fermions. In the process of propagation a particle can form 
bound states with other fermions according to the following mechanism. Some pair of fermions decays in components 
with momentums — k, —uj and k, uj with amplitude iA + . Second particle of the decayed pair is in state of the additional 
particle (k, uj) and it is identical to the additional particle. The second particle propagates through a system further. 
First particle of decayed pair forms bound state with the initial additional particle with amplitude — iA. Anew formed 
pair replenishes the condensate of pairs in a system. Thus, the dressed propagator Gs takes into account interaction 
of a particle, initially described by free propagator Go, with fluctuations of pairing. Intensity of the interaction is the 
amplitudes — iA and iA + . 

Starting from the aforesaid, we can write the mass operator for such process (FigH} as 

AA+ 

-iE = -iAiG (-k,-uj)iA + £ = (35) 



10 



-iZ = 




Figure 4: The diagram for the mass operator E describing interaction of a fermion with fluctuations of pairing. 

This mass operator has been proposed in 0, [2l[, however an existence of the amplitudes A, A + and the equation of 
self-consistency were postulated (as the anomalous averages). From Dyson equation (|A6[) we can obtain the dressed 
one-particle propagator: 

n f, s 1 u + e A s , B s 



w-e-S oj 2 -E 2 u-E lo + E 
G s (k, t) = -i6 t A s e- lEt + i0- t B s e iEt (36) 



where 



E(k) = \Je 2 (k) + A 2 (37) 



is dispersion law of dressed particles (quasi-particles). The amplitude A is named by gap, because minimal work for 
creation of one-particle excitations is 2A. In accordance with the definition of one-particle propagator we can write: 

G fk t) - I -^1^(^(0)1*0), t>0) 



where the system is placed in another ground state V&o- In the state interaction of particles with fluctuation of 
pairing (existence of condensate of pairs) is taken into account, moreover "Jo (A = 0) = $o, Gs(A = 0) = G . A 
propagator defines occupations number of quasi-particles n k by the following way: 

n k = (*o|C+C k |^ ) = -i lim G s (k, t) = B S , 1 - n k = (* |C k C+|* ) = t lim G s (k, t) = A s . (39) 

Hence, we can suppose, that 

C+|*o) = V^I*o, l k ), G k |* ) = v^l^o, l k >, (40) 

where \ ^o, 1^) and |\l/0) l k ) are states with one added particle and one removed particle with momentum k accordingly. 
Hamiltonian of a system of free quasi-particle has a view: 

k 

Chemical potential of a quasi-particles' system equals to zero. Therefore the grand potential fl coincides with 
Helmholtz free energy in a superconductive state. The states l^o), I^Ojl^.) and |^ojl k ) ar e eigenvectors of the 
Hamiltonian (1411): 



i/ |*o) =n |* ), #o|*o,l|;} = (fio + E(A0)l*o,l£}, ffo|*o,l k ) = (Oo + ^(fc))|*o,l k ). (42) 
Then, using the definition (|38|) . we can find: 

G s (k,i>0) - -»<*o |C k (*)Cj (0) | * > = -<<* |e i2ot Ci c e- i2ot Cj|*o) 
= l^|e <n ° t V^'e- <(nD+BW) *V^|*o 1 l k ) 

= -iAse""^*,,, 1P|* 0) iP ) = -iA s e~ iE \ (43) 
that coincides with (131)1) . For t < the proof is analogous. 
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Dyson equation can be represented in other form. Let's use the definition (ui — e)Gq = 1. On the other hand we 
have Go = Gs/(1 + Gs£). Moreover, let's introduce the notations 

-G S £ = AF+, -G S £ = A + F. (44) 

Then, we can obtain the set of equations: 

(w-#))G s (k, w ) + A_F + (k, w ) = l (45) 
(u + e{k))F+(k,u) + A+G s (k,w) = 0. (46) 

These equations are Gor'kov equations in momentum representation. However, unlike phenomenological approach 
(where existence of the anomalous propagator F and the equation for order parameter are postulated) these equations 
are obtained by microscopic way with help of the procedure of uncoupling of correlations. From the equations (I45I46[) 
we can find, that 

-A+ -A 
F+ = ^^2> F = ^JF2- (47) 

The anomalous propagators describe creation of two fcrmions from the condensate of pairs - F + , formation a pair by 
two particles with leaving to the condensate - F. Moreover, F and F + are the infinity sum of the serial processes of 
creation and annihilation of pairs described by amplitudes A and A + . Mathematically this is expressed in the fact 
that 



^(M) 



A+ f i<*o|C^ kp (t)Cj a (0)|*o>, t > 
Va+A \ ^o|G+ j(0)G+ k /3 (i)|*o), t < 
A+ ,—— .„ A+ 



g a p Uo t -^=^MB~ s e- iEt + i6- t -j==jMB~ s e iEt ) , (48) 

i(*o\C ka (t)C_ k p (0)\*o), *>0 
^A+A I i(*o|C_ M (0)C kiQ (t)|*o>, t<0 

A ?T~ -iEt , SD A n-5-AEt\ ( 1 



<? a/3 ( i6 t -^=yUsBse- iEt + t0- t ^==y/AjBse lEt ) , W = JJ). (49) 



Let's prove the formulas (|48I49|) . For this let's consider the expression: 

<*o|G+ k ^)G+ Q (0)|* o ) = (^ O \e iSot C+ kt0 e- iS ^C+J^ o ) 
= (*o, lV |e in °*V^e- i f ao+E( : fc »*V^|*o, l£ J 



= ^A^e-^^o, l h _ k J*o, l P k J = ^B~ s e-^\ (50) 

that corresponds to f|48|) ■ In the last equality the fact has been used, that the states, created by addition of a particle 
to a state k,a or removing of a particle from a state — k, j3 at a ^ /3, are identical in the course of existence of the 
pair condensate. The rest cases is proved analogously. It is not difficult to see, that F(A = 0) = F + (A = 0) = 0, 
because Ao(k)Bo(k) = 0. From (|48j) we can see, that the existence of nonzero anomalous propagators F and F + 
means breakdown of global gauge symmetry in a system, that is number of particles is not conserved in the course of 
existence of a pair condensate. Hence the states $o and have different symmetries: 

($o|*o>=0. (51) 

However distribution function over iV has a maximum at the average number of particles (N) determined by the 
expression: 

<*> _ 2 5> k _ -* / to O.fc t) - 2 / » N - 2 J (52) 



Now let's return to the left part of the expression (|25[) for T-matrix. After considering of particles' interaction with 
fluctuations of pairing we have to substitute dressed propagators Gs instead free propagators Go in the formula (|25l) . 
It is not difficult to verify, that T hasn't poles at any a and A. This means, that the problem of instability of a system 
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is removed, and we can use dressed propagator for further calculations confidently. Moreover, the absence of poles 
means the absence of bound states, because we have taken into account them in the specter of quasi-particles E(k) 

(37). 

It is necessary to note, that the mass operator (|35[) and Gor'kov equation (|45I46[) haven't parameter of interaction 
between particles. This means, that the amplitudes of pairing A and A + exists regardless of interaction between 
particles and its type. However, as we shall see below, the interaction determines the average value of the amplitudes, 
that is observed in experiment. This average value is not zero in the case of attraction between particle only. 



IV. GROUND STATE ENERGY. 



A. Summary kinetic energy of particles. 



In order to calculate ground state energy it is necessary to know kinetic energy of particles and energy of their 
interaction. The operator of kinetic energy of all particles of a system has a form: 



W = £ v F (k - k F )C+ a C k<a = 2 £ e{k)C+ C^ a . 



(53) 



k,a k 
Then the corresponding average value is 

—G(\<L,uj)e~ lut e{k) 
2tt 



/oo 
B{e)ede. 



(54) 



k _____ 

Here v F = ^ is density of states on Fermy surface. Since the interaction Vj_|j c _| c (|24p exists in the layer —uid < 
e(k) < ujd only, we can suppose that 



G 



G ; |e(fc)| > lo d 
G s ; \e(k)\ < uj d 



A(k) 



A (k); \e(k)\ > w D 
A s (k); \e(k)\ < uj d 



B(k) 



B (k); \e(k)\ >u D 
B s (k); \e(k)\ <uj d 



Then we can separate a normal part and a superconductive part of the kinetic energy: 

B a eds + Vv F j B s ede + Vv F J B ede 

WD fUD 



W n + Vv F 



= W n + V^- [ Lj'ij - oj D Juj 2 n + A 2 + A 2 arcsinh^ 



Bsede — Vv F 



B ede = W n — V 



LO D 



v F r° (e 1 e 2 



2 J- UD V E lei 



We can see, that the pairing leads to a loss in the kinetic energy. 



B. Vacuum amplitude. 



(55) 



(56) 



Since we took account of interaction of particles with fluctuations of pairing and we discovered that the ground 
state of a system |\&o) has other symmetry as compared with the initial state |$o)i hence the vacuum amplitude of a 
system can be written as: 



R(t) = (vf |C/(i-to)|*o) to =o-(*o|C/(i)|*o)e Wot 



E 

n=0 



dh... dt n ^ \T{H I (t 1 )...H I (t n )}\y ), 



(57) 
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where H](t) — e +ltH °V e~~ ltH ° is operator of particles' interaction in an interaction representation, Wo - ground state 
energy of a system without interaction. The averaging ()o = (^ol • • • l^o) is realized by ensemble of noninteracting 
quasi particles. The Hamiltonian of such system has the following form: 

fl„ + F = ^^#)C+ a C k|a + i Yl E V klm nC+ p C+ a C m ^Cn,s, (58) 
a k «,/3.7.5 k,l,m,n 

where momentum is conserved k + 1 = m + n and spin is conserved a + /3 = 7 + 5, V is volume of a system. The 
sequence order of indexes of matrix elements and of creation and annihilation operators is important. 

The expressions for several first orders (n — 0,1,2...) in the expansion of vacuum amplitude are (let us suppose 
i 2 > ti for definiteness): 

R (t) = (*q|*o> = 1 

R ^ = TiTf[ dtl E E (-^Vkimn) (*o|C+ (h)C+ (ti)C m , 7 (ti)Cn,i(ti)|*o> 
1 Jo ■ — kin... ' 

",/3,7,<5klm,n V 7 «'^',7'^'k',l',m',n' V 

x (*o|q^2)C+ Q ,(i 2 )^^^ (59) 

where C + (i) = e +4 * Ho C + e~ ltHo , C(i) = e +Itff °Ce" ,lffo are operators of creation and annihilation in Heisenberg 
representation, that coincides with interaction representation for an ensemble of noninteracting particles. 

In Appendix [B] we propose the method of uncoupling of correlations for approximate calculation of a vacuum 
amplitude R(t). As an example Hartree-Fock normal processes have been considered there. We shall generalize this 
method for anomalous processes here. In our case the particles interact by the potential (pM|) Vj 1 k _j c . Hence, the 
vacuum amplitude has a form: 



R{t) = 1 + W / ^EE (-H-l,k,-k) (*o|C^ lfl (ti)Cf a (ti)C kia (t 1 )C_ ki/J (t 1 )|* ) 

4^jf*jf *EE ("h-u.-k) EE (-K-r*-.V 



k,i x ' «'./9' k ',r 

x<*o|C+, ^2)^^*2)^ (60) 
We can uncouple correlations by the following way taking into account anticommutation of the operators C and C + : 



«./3 k,l V 7 

4<-^jf *jf *££(-5n-u-k) E E ( :>r .*■<) 

" , ° - 70 a,/9 k,l V 7 «'./3'k',l' V 7 

+ ... = l + i? 1 + ii? 2 + ... = cxp(i? 1 ) (61) 

As a result of the uncoupling of correlations we can see, that anomalous processes CC and C + C + give contribution 
to the vacuum amplitude of a system with the interaction (|24[) only. Another combinations of the uncoupling with an 
obtaining of normal propagator C + C don't conserve a momentum. Such representation of the vacuum amplitude by 
uncoupled correlations is analogous to Fock approximation for normal processes, and it means a neglect of dynamic 
correlation between pairs. Then R(t) can be written as 

bxR(t) = R 1 (t) = i r^^rW-ly^i.^ 

° «,/3 k,l V 7 
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where a summation over spin variables gave the multiplier 2. By analogy with Fig |12l in Appendix [Bl the process 
of uncoupling of correlations we can represent by graphically in Fig|5j where in the diagram of a scattering of two 
fermions from the states k and — k into the state 1 and —1 as a result of the interaction Vj j ^ _j c , we connect the 
lines corresponding to oppositely directed momentums. As a result we have two anomalous propagators: F(k,uj) 
and F + (1,oj). The obtained diagram means the following. Two fermions with opposite momentums and opposite 
spins appear from a condensate of pairs, interact with each other by the potential Vj _j ^ ^ and go back into the 
condensate of pairs. 




Figure 5: The uncoupling of correlations in a vacuum amplitude for the process of scattering of two fermions by the potential of 
interaction Vj _j ^ ^ from the initial state k, a and — k, /3 to the final states 1, a and — 1, (3. As a result we have the anomalous 
transition amplitude " vacuum- vacuum" . 

Using the formulas (|4"5|) and (|49p we can rewrite R(t) in the following form: 

\nR{t) = ^W^f) t y/A{e)B(e)de / y/A(e)B(e)de = -iXV ( ^) A 2 arcsinh (^J t. (63) 

— LJ £) — U £) 

This formula was obtained from such arguments. Since the area of action of the potential is limited by the layer 2ujd 
in a neighborhood of Fermy surface, then the amplitudes A and A + is not equal to zero in this area only. As it was 
pointed before, presented manner of uncoupling of correlations is analogous to Fock exchange interaction for normal 
processes. In Hartree-Fock approximation a decay of quasi-particles is absent jl|. This means, that the amplitude of 
pairing is real: A = A + in a momentum space. 

For calculation of the contribution of interaction to internal energy it is necessary to use the theorem ©, which 
connect a vacuum amplitude with a ground state energy: 

n A = I ii„fl( ( ) = -Ayj / |%f+ M £ / 

-"(t)v:^/:a— m 2 ^- (£). 

Since A < 0, then the interaction tries to reconstruct a system so, that the gap A is as much as possible. We can 
use the theorem because interaction of particle with pairing fluctuation, changing a symmetry of a system: $o — >■ \l/o, 
Go — > Gs, was considered before switching of the interaction Vj _j ^ ^. Switching of the interaction _j ^ jj, 
transferring the state "Fo in some other fy' , doesn't change symmetry of a system: (^ol^o) ^ This means that the 
adiabatic hypothesis is correct in the case of phase transition even. 
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C. Internal energy. 



On the basis of above obtained results we can write expression for internal energy of a system (at the temperature 
T = internal energy coincides with free energy): 

n = (W) + n x = -2ilna_J2^G(^u;)e-^e(k)-^Yl J ^ F+ ( k ' w )E / ^ k F(k,u). (65) 

k k k 

As it was pointed before, the gap is real A = A + , then F(k,ui) — F + (k, to). We can see, that the energy depends on 
the unknown amplitude of pairing A, which corresponds to Bethe-Salpeter amplitude r\ in the two-particle problem. 
The amplitude A is determined by dynamics of all particles of a system, and its observed value is a result of averaging 
over a system. The procedure of averaging means mathematically, that the observer value of A minimizes the internal 
energy: 

— =0^H)A=-£ J -w k F(^), (66) 

k 

that coincides with (|17[) . It means that the self- consistency equation for order parameter in Nambu-Gor'kov formalism 
is an extremal of the obtained free energy functional (|65]), and the order parameter is the averaged Bethe-Salpeter 
amplitude over a system. 

The functional (|65p can be written in an explicit form in quadratures: 

2 J_ UD \E \e\J 2 y J_ UD 2E J_ UD 2E V ' 

where VL n is the energy of a normal phase, g = A^ is the effective interaction constant. Value of the energy on the 
extremal (155)1 is 

fimin = V n + -UJ D ^ 2 D + A 2 ). (68) 

Thus, we solved the basic problem of statistical mechanics (at zero temperature): the calculation of a partition 
function (free energy) and, in particular, of a vacuum amplitude in a system of interacting particles for case, when the 
interaction causes a phase transition, that is symmetry of a system changes at a switching of the interaction. Unlike 
other methods, this result was obtained from first principles without introducing any artificial parameters of type of 
order parameter, but starting from parameters of the Hamiltonian only. 



V. NONZERO TEMPERATURES. 

A. Normal and anomalous propagators. 

In the sections IIIII and IIVI we have described microscopic mechanism of formation of long-range order in a system 
at zero temperature. In this section we shall formulate the approach for case of nonzero temperatures. Let we 
have a system from N noninteracting fermions being in volume V at temperature T. Then we must use Matsubara 
propagators, where time t is complex: t — > —ir, r £ [0,/3]. In ideal Fermy gas propagation of a particle with 
momentum k, energy e ~ Wf(|k| — hp) and spin a is described by the free propagator: 



G (k, r = T 2 - n) = 



-iSpjpoCk^C+^n)}, r>0 
iSpjftC^TOC^te)}, r<0 
-ie r {g+ A e-^ T + goB e^ T ) + A Q e-^ T + g+B e^), 

n=+oo „p 

G(k,r) = ~ G(k, W „)e-^ T , G(k,w n ) - - / G(k,r)e— T dr 

" n=-oo 

G (k,w n ) = - 7T\= l T — \2 2 =l ' n +l ' TT7' ( 69 ) 

iuj n - e(k) (iLo n ) 2 - e 2 iuj n - e iuj n + e 
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where 

,1 1 (2n+l)7r 

50 = ^TI' % = iwJ+T 1 " n = (70) 

Cjj. ^ (r) and G^ (r) are operators of creation and annihilation in Heisenberg representation, po is density matrix of 
noninteracting particles: 

in-Ho+jm) f "-E k , g ^, g c£, g + /^ l ( "-Sk^W^^ l ^ 

Po = cxp <^ S = exp ^ i V = exp <^ '- S , (71) 

where e(fc) = j— — p ^ vp{k — kp) is kinetic energy of particles counted off from Fermy surface. 

Now let an attracting force acts between particles. The force is described by the matrix element of interaction (|24l) . 
In this case the instability of a system appears again with regard to a pairing (as in the section [II]) and T-matrix has 
a pole structure: 

1 (U,U) — — r- 2 — — — , Ic= — Wd I -^TT-: , (72) 

l + A^^fln^ 2 - uiKfT-Tc 7t \ 2|A|j/f/ 

where i/p = ^^r 1 is density of states on Fermy surface. We can see that bound states exist in a system while 
temperature is not higher than critical temperature Tq - at higher temperatures particles have large kinetic energy, 
so that an attraction between them leads to a scattering only. 

Generalization of the two particle problem on the multiparticle case is done analogously to the section Mil In the 
case of nonzero temperature the mass operator has a form: 

-AA+ 

- E(k >Wn ) = (-A)iG+(-k,u; B )(-A+) = (73) 

lUJ n -f- £(K) 

Then it follows from Dyson equation, that a dressed propagator has a view: 

11 i 

iE^>G s (k,w„) 



G G s " ' iu) n - e(k) - E(k,w„) 

iw n + e .A s B s , 
»r [— r + «r j— r- (74) 



(iw„) 2 - £ 2 (fc) ioj n -\e\ iu) n + \e\ 
It can be written with help of a total definition of Green function in (k, i)-space: 

| l Sp{g C+ ff (r 1 )C kiff (r 2 )}, r<0 

= -iB T {g^A 8 e- BT + 9 S B s e Er ) + il T ( 9s 4 s e- ET + g + s B s e ET ), (75) 

where 

^ = ^w~rv g s = ^ftt- (76) 

go is the density matrix of noninteracting quasi-particles: 



*n-Ek,^(*)Ck,^ 



g = cxpj ■ — ■ (77) 

It should be noted, that the state described by the density matrix go has another symmetry in comparison with the 
initial state pq. Occupation numbers n(k) is determined by the following manner: 

n(k) = -i lim_ G(k, r) = g~ A + g + B, Urn n(k) = B(k). (78) 
Let's introduce the designations: 

-GE = AF+, -G + S+ = A+F.. (79) 
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Then Dyson equation can be rewritten in a form of Gor'kov equations: 



(iw n - e)G + AF+ = i 
{iu n + e)F + + GA = 0. 

The expressions for anomalous propagators follow from Gor'kov equations: 

-iA+ 



{iw n y-E*{ky 



F(k,u„) = (F+(k,w„))- 



iA 



(80) 
(81) 



(82) 



These expressions are analogous to the expressions (|47|) . We can write the anomalous propagators in (k,t)-space and 
in a form of a vacuum average of creation and annihilation operators: 



^(k,r) 



A+ [ »Sp{goC+ kff (T 2 )C+ a (n)}, r>0 



A+ 



^(k,T) 



'Va+a 

that is analogous to the expressions (14"51) and (|4l?)) . 



A f -<Sp{gbC kia (7a)C_ kij9 (r a )} J r > 
I -*Sp{aC_ k ^(ri)C k a (r 2 )}, r < 

igap-^y/A^ [- (g+e- ET - g~ s e ET ) T + (g+e BT - ^e"^) 0_ T ] 



(83) 



(84) 



B. Kinetic energy and entropy. 

In order to calculate a free energy it is necessary to know kinetic energy of particles of a system, energy of their 
interaction and entropy of collective excitations. Average kinetic energy of all particles of a system is 



n— +oo 



(W) = -2i5^G(k,T->(r)e(k) = — ^ lim J2 G(k, w n ) e -^"*e(fc) 
k T ^ k ™=-°° 

/OO 
{g-A + g+B)ede. 

Since the interaction ^_ljc_k exists only in the layer — tj£> < < cjd, then we can suppose that 



(85) 



.9o ; H fc )l > u D 

9s~\ I £ ( fc )l < U D _ 



g + = 



9o\ |e(fc)| > 
5 +; |e(fc)| < w D 



(86) 



Hence, one may write expression for kinetic energy separating normal and superconductive parts by analogy (|56[) : 



(W) = »'„ - V^f f° (g~ - „+) 



- r ° ' - « ids 



[9s -9s, 



E 



(87) 



In the limit of low temperatures /3 — » oo this expression reduces to the expression (1561) . If we suppose that A = 0, 

then we shall have W = W n . 

At temperature T / a gas of collective excitation exists - boholons with the spectrum E = ^e 2 {k) + A 2 . Since 
boholons are product of decay of Cooper pairs on fermions, hence occupation numbers of states by boholons are 



M fc ) = ^i = U 1 - tanh (f) 
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Then entropy of a system is 

S = -2 [/(*) l n /(*) + (!"/(*)) Ml "/(*))] 
k 

= S -2V-f [fskxfs + (l-fs)Hl-fs)}de + 2V-f [/ ln/ + (1 - /„) ln(l - /o)] de, (89) 

Here we separated the normal part again, where /o = (e' 3 ' 6 + so that S = S n at A = 0. The multiplier "2" 

appeared as result of summation over spin states. 



C. Vacuum amplitude. 

In the previous subsections we considered the interaction of particles with fluctuations of pairing, and we found, 
that the state of a system described by the density matrix go has another symmetry in comparison with the initial 
state pq. Hence vacuum amplitude of a system can be written in a form: 

R((3) = (U(0)) o = Sp (qoU(0)) = —T / dn... dr n Sp (g T {fi^n) . . . ff/(r„)}) , (90) 

where H](t) = e + TH °Ve~ THa is the interaction operator of particles in interaction representation. The averaging 
Sp (j>oU(/3)j is made over ensemble of noninteracting quasi-particles. The potential Vj j k k acts between particles 
Hence we can write (by analogy with (|60[1 ) expended expression for vacuum amplitude: 



R W = 1 + 77T7 /" dTl EE (-^l.-lk,-k) Sp{&C+ 7 (T 1 )C+ a (r 1 )C' kia (r 1 )C_ ki7 (r 1 )} 

^^f^f^EE E E ( 2 «i.i*.v 

J0 J0 a,7 kl V 7 «',7'k' 1' 



+..., (91) 

where we took into account that Sp{go} = 1- Then, by analogy with the equation (|61|) . we can uncouple correlations 
by the following way taking into account anticommutation of operators C and C + : 



m « 1 + 7 (-i) 2 ^ / d n £ E f-H-l,k,-k) sp {^oC+JtOC+j (n)} s P {s,c_ k , 7 (n)C k Q ( n )} 

k.l V J 

a,<y k,l a'^'k',1' 

xSp {ebC+ Q ,(r 2 )C+, >y (7a)} Sp {ftC_ k . r (T 2 )C k . a ,(r 2 )} Sp {^C+JrOC+^n)} Sp {?oC_ k , 7 (ri)C7 k Jn)} 
+ ... = l + Ri + -Rf + ...=exp(R 1 ) (92) 

Let's take into account that our approximation is analogous to Fock approximation for normal processes. Then we 
can suppose A = A + , hence F = —F + . Then R(t) can be written as 



7 J Q dTl E E {^\ v l-lX-kj iF °<i( l > Tl _ n)iF aP (k, n - tx) 



a 'f k.l 



V E (Vl, lk k) *"(!. - 0-)F(k, r -> 0-)/3 



v ^ \2 

k.l 
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^ n— -\-og n— 4-oo 

|^ 71 — — DG |^ n — — oo 



■^(f)7::-(f)^/:>-(f)^ 



This equation is represented graphically as well as in Fig[6j In order to calculate a contribution of interaction in free 
energy we can use the formula (fTU)) : 



. , n— +00 1 +00 

k n— — 00 ]^ n— — 00 

-nf)7:-(f)^/:-(f)^- 



In the limit of low temperatures — > 00 this equation transforms to the equation (|64[) . If we suppose that A = 0, 
then we shall have fix = 0. 



D. Free energy. 

Starting from the above found results we can write the expression a for free energy of a system: 

r> ■ n=+oo 

n = {w)-~s + n x = -- um V V G^^y-^eik) 

' 71— — 00 

f~. , 1 n— +00 -j n— +00 

+7Ei/w in /w + ( 1 -/w) ln ( 1 -/w)]-^E w '« E ^(k,^)E^-« E F ( k > w «) ( 95 ) 

n— — 00 ]^ n— — 00 

We can see, that the energy depends on the unknown amplitude of pairing A, which corresponds to Bethe-Salpeter 
amplitude r/ in two-particle problem. The observed value of A must minimize the free energy: 

g = => (-i)A = ^ E E ^ k ' 06) 

n— — 00 

that coincides with (p~7j) . As in the previous section the equation of a self-consistence for order parameter in Nambu- 
Gor'kov formalism is the extremal of the functional of free energy (|95p. and the order parameter is averaged Bethe- 
Salpeter amplitude over a system. 

The functional (|95p can be written in an explicit form in quadratures: 



, . /3 e \ e 2 , //3£\ e 2 
tanh ^ — - tanh - — — 
2 e \ 2 J E 



de 



+yy| D [fs In /s + (1 - fs) ln(l - /s) - /o In /o - (1 - /o) In(l - /„)] de 

where £l n is the energy of a normal phase, 5 — X^- is the effective interaction constant, V is volume of a system. 5 
can be expressed via critical temperature /3c with help of the equation A (/3c) = as following: 

1 = - 5 /"" D tanh ( ^) -Lde. (98) 



— 



2e 



If to suppose A = 0, then we shall have 51 = tt„ . The equilibrium value of A is determined by balance of kinetic 
energy, entropy and energy of interaction, that corresponds to a minimum of the free energy. 

Let's consider a low-temperature limit of the free energy (|M)l : A/3 > 1 at T -> 0. This means, that the value 
A — Aq can be parameter of expansion, where Aq = A(T = 0) is equilibrium value of the gap (amplitude of pairing) at 
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zero temperature. Due a rapid convergence of integration elements in (|96|) . the limits of integration can be — oo, +00. 
Then a low-temperature expansion has a form: 

Q = n n + V (a (T) + b (T)A + d A 2 ) , (99) 

where coefficients of the expansion are 



a (T) = _^ + ^-^v^^e-^+^( ff + l)Ag-^( ff +l) ) /ftrA8re-^ 
ho CO = -^2( 5 + l)A + ^( 5 + l)v^A^r e -^, d = ^2( 5 + l) (100) 

Since in the expression ^rf the multipliers 5 + 1 is cancelled and, as a rule, |<?| <§C 1, we can suppose that g + 1 « 1. 
Then expressions (|100[) are simplified: 



ao(T) = ^ + 2, F T 2 -^V^^e-^+-^ v /8^T e ^ 
6o(T) = -^2A + ^v/8^A T e -^, d a = v F (101) 



n = ft n + V I a(T)A 2 + j-bA 4 + ^dA 6 ) , (102) 



Let's consider a high-temperature limit of the free energy: Aftc <C 1 at T — > Tp. Expansion in powers of A gives 

i&A^ + i 
2 3 

where the coefficients of the expansion are 

v F T - T c 



a(T) 



2 T r 



b 2 8^ 2 T 2 [W6) 
rf ^/52.31 C (5) +4 ^ 1 



4!T 4 ' 



This expansion has a form of Landau expansion of free energy in powers of order parameter. 

From the all considered above we can see, that averaged over a system Bethe-Salpeter amplitude r\ and n* - the 
amplitude of pairing A and A + have the properties, which are analogous to the properties of a order parameter. Thus, 
we have solved the basic problem of statistical mechanics: the calculation of a partition function (free energy) and, in 
particular, of a vacuum amplitude in a system of interacting particles for case, when the interaction causes a phase 
transition, that is symmetry of a system changes at a switching of the interaction. 



VI. THE PAIRING WITH NONZERO MOMENTUM OF CENTER OF MASS OF A PAIR. 

Let fermions with momentums k + 9 and — k + 9 pair up, so that the momentum of a center of mass of a pair is 
q. The free propagators, corresponding to these states are 

G 



2 > c-e(k+§ 



uj — £4 



G (-k+3_a,) = j _ = — -. (104) 

2 ' -w-e(-k+SP -u-e- 



We can suppose that 



e+ = TT-fk+5 N ) 2 -uwe+-5-. e_ = J- f-k + ^) 2 - a rj e - (105) 
+ 2m V 2 / P 2m ' 2m V 2 / P 2m V y 

The corresponding mass operator describing an interaction of particles with fluctuations of pairing is: 

(-i)E, = -iAiG f-k + ^, -a;) ^A+ = (-i) AA+ . (106) 
V 2 / cj + e_ 
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Then can find the dressed propagator from Dyson equation: 

1 _ w + e 
Go 1 - S g ~ (u>-E+)(u - EJ) 

where the specters of quasi-particles are 



£ + f £ + +£ -\ , I A 19 kf l 



iL-ayz-^^j + |A|- S -V? + W. (108, 

We can see, that if to assume q = 0, then the specter (|108[) turn into the usual specter of boholons: E = ±\/e 2 + |A| 2 . 
The critical momentum g cr exists when a minimum of the specter( |108|) touches Fermy surface. Then for excitation 
of a system it is necessary infinitely small energy. This means, that superfluidity of Fermy gas is absent. The critical 
momentum is: 

E+(q = q cr ,e = 0) = E_(q = q cr ,e = 0) = q cr = — |A|. (109) 



Similar pairing can take place in high-temperature superconductors (cuprates) with a mirror nesting e + 3 
e(— k+3j f regions of Fermy contour (2(| [27j • 



VII. FREE ENERGY IN A CASE OF SLOW SPATIAL INHOMOGENEITY. 



In the previous sections we supposed, that amplitudes of pairing A and A + don't depend on spatial coordinates. 
This takes place in interminable, homogeneous, isotropic and isolated from external fields superconductor. However 
in a total case these conditions are not realized. For example, in a sufficiently strong magnetic field the inclusions 
of normal phase can exist in volume of a superconductor. Another example - a contact of a superconductor and a 
normal metal. In this case the order parameter is suppressed in a boundary layer of a superconductor, however it 
appears in boundary layer of a normal metal. 

Let's consider some region of a superconductor, where a distribution of A is inhomogeneous and the amplitude 
of pairing can be both smaller and larger than its equilibrium value Ao - Fig[SJ When a quasi-particle propagates 
along the axis Ox its energy is constant E — \J e(k) 2 + | A| 2 , but its momentum changes. If a pair moves into region, 
where |A| < |Ao|, then forces appear tearing the pair. Each element of the pair gets some increment of momentum q: 
(k, — k) — > (k + q, k — q). If the pair moves into region, where |A| > |A |, then a forces appear increasing a bound 
energy of the pair. Each element of the pair gets some increment of momentum q too. Since the amplitude of pairing 
is function of coordinates A(r), moreover we suppose that order parameter is real A = A + in momentum space and 
it has an identical dimension in q-space and in r-space, then we can write the Fourier-transformations: 

A(r) 
A+(r) 
A(g) 

The mass operator for above-mentioned process is shown in Fig|7] In analytical representation it has a view: 

- S,(k,a; n ) = (-A,)*G+(-k- q,u, n )(-A+) = -r , (HI) 

where the free propagator Go is 
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Figure 6: The pairing of fermions in a spatially inhomogeneous system. Ao is the equilibrium value of a gap in a homogeneous 
system. 



iG n (-k-q,co n 



Figure 7: The diagram for the mass operator £ describing an interaction of a fermion with fluctuations of pairing in a spatially 
inhomogeneous system. 



Then from Dyson equation we can obtain the dressed propagator: 

— - — »t =4 r? c -v l ^ + £ i 

Go - g s ~ Gs ~ Vn) 2 - 



where E q is the specter of quasi-particles in a nonhomogeneous system: 



El = el + \A q \\ 



kq 



e q = e(k + q) ss e(k) H , |k| ~ Hf- 

m 



(113) 



(114) 



Dyson equation can be represented in a form of set of Gor'kov equations, from where the expressions for anomalous 
propagators follow: 

(tw n - s q )G + A q F+ = i F+(k + q,cj n ) = { .J n )tl E , 

=> (H5) 
{iu n + e q )F+ + GA q = F(k + q5 Wn) = (F+(k + q> w „))+ = {iM lfi_ E , 

If to suppose 9 = 0, then we shall have the expressions (|80M82I) . 

Now let's suppose that A(r) changes very slowly on a coherence length l(T) which characterizes a size of Cooper 
pair. Then we can suppose e(k + q) = e(k) in the specter of quasi-particles, such that E q w yje 2 {k) + |A g | 2 . However 
we must keep e(k + q) in numerator of the expressions (|112l) and (111 3[) for G. Hence the normal propagator has a 
form: 

G(k + q, r) = -iB r (g+A(k + q) e - B * r + g q B(k + q)e^ r ) + i0_ T (g q A(k + o^e'^ + g+B(k + q)e^) , (116) 
where 



A(k + q)« ^(fcH-L^, fl(k + q)«S ff (fc)- ' k< " 



2£ , (J m 



The anomalous propagators are 



F+(k + q,r) = iflafl^- e 



F^(k + q,r) = i ga/} £±[-(g+ e - E ' r 



2E q m 



(g+e E ^-g-e-^)e. T ], 



(g+e^-g^e-^) 



(117) 



(118) 
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We can see, that in the approximation of slowness of changes of A(r) the anomalous propagators depend on q by 
means of A(q) only. 

Kinetic energy of a system is determined by the following way: 

(W) = -2i^ £ (k + q)G(k + q,T^0-)=2^ e (k + q)( 9 -^(k+q) + 5 +S(k + q)) 
k k 

1 (kq) 2 



= 2 J2 £ ( fc ) (9q A q + 9tB q ) + 2 ]T (<?" - g+) 
k k 

^ 2 (PE q 



E q m 2 



W n + V 



Tr V F 1 2 2 



U1 D 



tanh 
tanh 



2 
2 



i 

W\ 



— tanh 

— tanh 



V 2 



1 

~E« 



de 
de. 



(119) 



We can see, that the term, which is proportional to q 2 , is added to the kinetic energy (|57|) (with the replacement 
E — >• Eg = ^/e 2 (k) + |A(<?)| 2 ). In the approximation of slowness of changes the expressions for entropy and vacuum 
amplitude coincide with the expressions (j8"9")l and (IM1) accordingly, however it should be written A(q) instead of 
A = const. Then we can write the free energy: 



n( g ) = n n ( g ) + n(A q ) + V~v 2 F q 2 



tanh 



2 



— tanh 



0E q 



1 



de, 



(120) 



where fi(A g ) coincides with the expression (|96p . where the replacement A — > A(q) was done. 
Expanding the free energy (|120[) in powers of A we can obtain the expression: 



fife) = n n (q) + V ( a(T)A q + -bA q + jq 2 A 2 



(121) 



where the coefficient 7 is 



7 



7C(3)4 
2 24tt 2 T 2 



(122) 



where Zq is a coherence length at T = (Pippard length). The expansion (|121j) has a form of Landau expansion of 
free energy in powers of order parameter at the condition qlg <C 1. We can see, that a spatial inhomogeneity increases 
the free energy of a superconductor. Hence in most cases we can be restricted by the term ~ q 2 in the expansion, 
because more fast changes of A increase the free energy essentially. 

The full free energy of a system in a spatial inhomogeneous case can be obtained by summation of the expression 
(|120p over all possible q: 



y r Y 2 f 

fi = fi„ + =^n + ^3 J tts(q)d 3 q = fi„ + ^3 J 



a(T)A 2 ^6A 4 + 79 2 A^K', 7 



(123) 



Let's pass from momentum space to real space using the expressions (| 1 10p : 



A q Agd 3 q = / A 



A(r)e-^ r d 3 



d 3 q=y J A(r) J A{q)e~ lC ^d 3 



d 3 r 



(2tt) 3 
V 2 



J A(r)A+(r)d 3 r,(124) 



qA 9 qA 9 A = J qA 9 i j {-i)^A{v)d 
qA(q)e- ic Fd 3 q 



V 



d a r s ,3 ( 27 3 



d 3 g 

d . , s 

7T A r 
or 



^-A+(r)d 3 r. 
or 



(125) 



For the term A 4 and terms with more high powers the situation is more difficult. This is because a square of a Fourier 
transform is not equal to a Fourier transform of a square: (y J A(r)e _l( l r d 3 r) 7^ y J A 2 (r)e~ lC l r d 3 r. Apparently 
this fact results to some nonlocality of a superconductor's state in zero magnetic field described in [281 ] . where a value 
of gap in a point is determined by a distribution of gap in some neighborhood: A(r) = fdv'Q(r, r')A'(r). However 
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in first approximation this correlation can be neglected and we can write the expansion of free energy in powers of 
AA + in real space: 



ft = ft„ - 

This expansion coincides with Ginzburg-Landau expansion in zero magnetic field 



a(T)|A(r)| 2 + -|A(r)| 4 + 7 



(126) 



VIII. FREE ENERGY OF A SUPERCONDUCTOR IN MAGNETIC FIELD. 



In this section we shall generalize the previous results for the case, when a superconductor is placed in a magnetic 
field H(r) = rotA(r). Our aim is to obtain the functional of free energy ft (A(r), ^A(r), A(r)), which is correct for 
an arbitrary value of the relation A(T)/T, for an arbitrary scale of a change of A(r) in comparison with a coherent 
length Z(T), for an arbitrary value of a magnetic penetration depth A(T) in comparison with a coherent length Iq 
(nonlocal electromagnetic response). Thus, the all three restriction on Ginzburg-Landau functional, described in 
section U are excluded. 

Let the microscopic magnetic field exists in a superconductor with a potential A and an intensity H: 

A(r) = J2 a (qK qr =>" H ( r ) = rotA(r) =^qx a(q)e iqr . (127) 

q q 

Then the magnetic field inducts a current: 

J(r) - ^rotH(r) = E x x = ~T E ~ ac ?) = E Hl)^ ■ ( 128 ) 

11 71 q n q q 

Energy of the magnetic field is 

W f = ±j |H(r)| 2 d 3 r=^E tf< - (^) 2 ) ■ ( 129 ) 

A magnetic field affects on a superconductor essentially. In the first place, a distribution of order parameter becomes 
inhomogeneous. As it was shown in the section IVII1 an inhomogeneity leads to same growth of momentum of each 
element of a pair: k > k + q, — k — > — k — q, moreover the order parameter depends on the momentum A = A(q). 
In the second place, an ordinary momentum must be replaced by a canonical momentum: k + q— >k + q— |a 9 , 
moreover the order parameter depends on the momentum A = A(q — - a g ). 




Figure 8: The diagram for the mass operator E describing the interaction of a charged fermion with fluctuations of pairing in 
a spatially inhomogeneous system situated in magnetic field with a potential a(q). 

The mass operator for a process of interaction of a fermion (with charge e) with a fluctuation of pairing is shown 
in Fig(S] In analytical representation this diagram has the form: 

- S (k + q - -a„u„J = -A (q - -a g ) iG^ (-k - q + -a g , w n J (- A + (q - -a q 

- -Mq-klA+h-S-.), ' ' (m 

ioj n + e (k + q - |a g ) 



25 



where the free propagator Go is 

Go 



e(k + q- | a,) 



:(k + q-fa,) (iw n f - e 2 (k + q - § a,)' 



Then from Dyson equation we can obtain the dressed propagator: 



e (k + q - f a,) 



11. . 2W r , 

Go = G; _Zl]9 ^ Gs = Z (^0 2 -i? 2 (k + q-fa g )' 
where E is the specter of quasi-particles in a inhomogeneous system situated in magnetic field: 

I 2 = ^qa 



(131) 



(132) 



E z [k + q - -a 9 J = e a (k + q - -a, J +\A(q- -a q 

t e \ , , k (q — -& ) 

e(k + q--a,W fc + V c ~ = eq.a, k 

V c / m 



(133) 



where we have introduced the notations £q,a and eq.a for convenience. Then Dyson equation can be represented in 
the form of Gor'kov equations. From these equations the expressions for anomalous propagators follow: 



(iw„ - eqa)G + A q a F+ = 



F+(k + q-fa g ,w„) = ( , tL jf_ q j 



qa 



iton + eq*)F+ + GA q , a = p(k + q _ f ^ ^ = (F+(k + q _ ^ ^ ))+ = _jAq^_ 



(134) 



If to suppose g = a = 0, then we shall have the expressions (|80H82p 
In the space (k, t) the normal propagator has the form: 



G(k + q- -a 9 ,r) 



-*fr ^.<a^q,ae- £ q aT + gq , a Bq,a e ^ ^ 
»-r (ffq,a^q,a e - Bc l ar +ff+i?q,a e ^ ar ) , 

(k + q 



where 

A (k + q- -£ 

V c 

The anomalous propagators are 



1 I x , 



2J5, 



q : a 



q,a 



1 / x £q,a 



E, 



qa 



(135) 



(136) 



--Br 



^(k + q- -a„,r) 



iq,a 



2J5, 



qa 



g+ a e- E W r 9^ a e E ^ T ) 9 T + (g^ a e E ^ r gq, a e~ E ^ T ) 6^137) 



where g q a and g q a are statistical multipliers: 



9 



(k + q 



( k+q -^) = e -^q.a + 1 - 



(138) 



We can see, that the normal and anomalous propagators have a complicated dependence on vector q, amplitudes of 
pairing A, A + and magnetic field a(q). 

Kinetic energy of a system is determined in the following way: 



(W) = -2i^e(k + q- -aJ G(k + q 
k 

= E^ a fi- eq ' a 



-a q ,T -> CT 



2 X! £( l> a (ffq,a^q,a + ffq.a-^q.aj 



k 



tanh 



v q,a 



£ q a 

k,(|e(fc)|<u D ,|k|~fcjr) 



/ £q,a , /3|eq,a| £q,a , /3-Eq 

to tanh — — ^ 



= W n + W S . 



(139) 
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We can see, that the kinetic energy depends on the vectors q and a(q) in a complicated way. If we suppose a = 
and q is small, then we shall obtain the expression (|119p . 

Free energy of a superconductor is the sum of the following terms: 



Q = n n + W s - -pS s + n x + W f ieid(a), 



(140) 



where Ws is the kinetic energy of fermions of a system in superconductive phase, Ss is the entropy of boholons, 
fl\ = —jj InR(P) is the energy corresponding to an interaction, Wf i e ld( a ) * s the energy of the magnetic field (11291) . 
The expressions for S and fl\ are obtained from the expressions (|89[) and (|94|) (without transition to integration over 
e only) by the replacement A — > Aq a , E — > -Eq.a, £ — > £q.a, / /q.a- Hence the free energy of a superconductor 
for given q is 



f}(q, a 9 ) 



fi„(q,a 9 )+ £ qa 

k,(|e(fc)|<wr>,|k|~fc F ) 

2 



£q,a , /3En,a 
tanh 



|£q,a| 2 £^q,a 2 







E 

k,(|s(fc)|<a> D ,|k|~fcjr) 
A 

V 1 . ^ . . 2E, 



E 

k,(|e(fc)l<"r>,|k|~fc F ) 



/^a ^ jq\a + (1 - /|,a) 1^(1 - /q\a) - /qa In /q,a - (1 - Jq,a) - /q, a ) 
Ana . 0E a . a ^ A, 



q- a tanh^ 

qa 



E 

k,(|e(fc)|<W£,,|k|~fc F ) 



2E, 



^ a tanh^ 



qa 



2 2 



where 



(qa,) ) = n„(q, a q ) + fi s (q, a,) + w f ie id(q> a 9 ), 



/q,a(fc) 



/q,a(*0 - fl £qa 



^q al 



(141) 



(142) 



are occupation numbers of states by boholons, moreover with help of the term /q a a normal part of entropy is 
separated, such that fis(A = 0) = 0. The full free energy is sum of the expression (|14ip over all possible q: 



Q = fi„ + { fi s(q, a g ) + w f ieid(q= a ?)} • 

q 



(143) 



Unlike Ginzburg- Landau functional the obtained functional of free energy (| 143[) is correct for an arbitrary value of 
the relation A(T)/T, for an arbitrary scale of a change of A(r) in comparison with a coherent length l(T), for an 
arbitrary value of a magnetic penetration depth A(T) in comparison with a coherent length Iq - it describes a nonlocal 
response to magnetic field. However the obtained expression is complicated for analyze. For its simplification let's 
suppose that A(r) changes in space slowly. Then it is necessary to expand the expression (| 141[) in degrees of q— |a g 

keeping terms which are proportional to the vector in second degree only. Then supposing £q,a ~ \Je 2 (k) 
we have: 
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(144) 



where g — A^y- is the effective interaction constant, and the gap is A = A (q — -a g ). 

Let's consider a high-temperature limit of the free energy: A(3 C <C 1 at T —> T c . This means, that the expression 
(I144p can be expended in degrees of Aq a : 



n 



n n + Vj2 (<*(T) Aq, a + ^A q , a + ±dA q , a + 7 (q 
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J2 Wa\ (qa 9 ) 5 
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Figure 9: The screening of the external magnetic field H ext , impressed in parallel to a superconductive cylinder, by inducted 
closed currents. The inducted currents are directed so as to compensate the external field. The resultant current goes around 
a lateral area of cylinder, but strives to zero inside of the cylinder. The area where the resultant current is not equal to zero is 
surface layer with thickness ~ A. The magnetic field penetrates in a superconductor in the depth ~ A too. 



where the coefficients ev(T),&, d are determined by the formulas (I103p . and the coefficient 7 is determined by the 
formula (|122p . The expansion (|145j) has a form of Ginzburg-Landau expansion of free energy in degrees of order 
parameter. Observed configuration of the order parameter A q ,a and the magnetic field a(q) minimizes the free 
energy: 



— = a(T)A q ,a + 6A q , a + dA q , a + 7 (q - -a, J A q . a = (146) 
AO e 2 

^ = =4- j(q) = 2e 7 qA q!a - 2 7 -A q!a a(q), (147) 

where j(q) is Fourier component of a current: 

c c 
j(q) = - — q x q x a(q) = -— (q(qa) - aq 2 ) . (148) 

These equations take more simple form in a transverse gauge q • a(q) = 0. This gauge gives a condition of closure 
of a current (special case of conservation of charge) : 

q ■ a(q) = => j(q) = --^a(q)g 2 => q ■ j(q) = & divJ(r) = 0. (149) 

47T 

The closed currents (|149j) screen a magnetic field in a superconductor - Fig|9] The currents is analogy to molecular 
currents of Ampere, their resulting gives rise to observed magnetic effects. In transverse gauge the functional of free 
energy has a form: 

fi = fi„ + V J2 (V)Aqa + ^A q .a + ^ A qa + 7 (q 2 + ^4) A qa) + ^ E <? a l ( 150 ) 
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The equations of extremals are 

a(T)Aq, a + 6A^ a + dA q>a + 7 (V + ^a^ A q , a = (151) 

j(q) = -2 7 ^A qia a(q). (152) 

From (|152p one can see, that the value Q = — 27^-Aq a is a Fourier transform of a kernel in the integral law of a 
magnetic response (Pippard law). The order parameter is function of q and a(q): 

AUT) ^)) - HQl (. - H T) + UQ) , (153, 

where smallness of 1/q in comparison with the coherent length l(T): ql(T) <C 1 is supposed, and we assumed that the 
coefficient d = for simplification. From the formula (|153[) one can see, that the kernel of the magnetic response Q 
is function of magnetic field. Hence the electrodynamics of a superconductor is nonlinear. If to suppose A = const 
at given temperature, then we shall obtain London equation: 

^ = - 2 4A 2 (T)a(q) , -^a (q) A 2 ( T) = (154) 

where A 2 (T) is the magnetic penetration depth in a superconductor. It is necessary to note, that A oc = _ 
This means that the exchange of mass, charge and concentration of superconductive electrons to corresponding values 
of Cooper pairs doesn't change the observed values. 

However it is necessary to note, that in a high-temperature limit a gap (and a kernel Q) depends on vector q and 
field a(q) strongly. From the formula (|153[) one can see, that the gap decreases at an increase of q, hence a magnetic 
penetration depth A increases. Moreover with a rise of temperature this dependence becomes stronger (at T = Tc 
we have A = 00). Besides at temperature T — Tc the critical magnetic field is zero He = 0. This means, that in 
the limit T — > Tc any magnetic field H(q) can not be considered as weak. Therefore it suppresses order parameter 
essentially and penetrates in a superconductor deeply (in macroscopic distant even). For example, the penetration of 
magnetic field along a core of Abrikosov vortex in a type II superconductor. Such structure exists in infinitely weak 
magnetic field at T — > Tc. 

For research of the nonlocal characteristics of the functional of free energy (|141I143[) let's consider a low-temperature 
limit A/3 ^> 1 at T — > 0. A value of gap is close to the value at zero temperature A(T) < Ao- Moreover, magnetic 
field is weak, such that it changes a value of gap lightly, that is the magnetic field is much smaller than critical field 
H <C He- Either as above, we assume that a change of a gap in space is slow. Starting from aforesaid and using the 
expansion (|99p we obtain the free energy: 

n = fi„ + V ("OCT) + &o(T)A q , a + rf Aq ja + 7 (q - ^a,)' A 2 ^) + ^ E (? a l ~ fa 8 *) 3 ) > ( 155 ) 

where coefficients cto(T), bo(T), do are determined by the formulas (|101|) . and coefficient 7 is determined by the formula 
(|122p . The observed configurations of order parameter A qa and magnetic field a(q) minimized free energy: 

80 / e \ 2 

— -0 => k (T) + 2d Aq, a + 7(q--a 9 J A q , a = (156) 

SO e 2 

T = j(q) = 2e 7 qA q . a - 2 7 -Aq ia a(q). (157) 

oa M c ^' 

In the transverse gauge q • a(q) = the functional of free energy and the equations for the extremals have a form: 

n = n n + Vj2 ("OCT) + 6o(T)A q , a + d Aq. a + 7 (q 2 + ^Ja 2 ) A q>a ) + £ E f a \ ( 158 ) 
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2 



b (T) + 2d A q , a + 2 7 ( q 2 + ? a 2 ) A qa = (159) 
j(q) = -2 7 ^-A q ^ a a(q). (160) 



29 



If in the equation (|160l) to assume A = const at given temperature, then we shall have London equation again: 

J(<0 = -27^(T)a( q) , -^«<q) - A»(D = ^^fp (161) 

The set of equations (| 159|160[) allows to generalize London equation. From the equation (|159j) we can find value of a 
gap: 

Aq a(T) = ~ b ° (r) 2 - (162) 

q,{U ' 2d + 2 7 (q 2 + #a2) V ' 



Then the equation for current has a form: 



e 



2 



b 2 (T) 



j(q) = -27-- U ' - „ 2 a(q) = Q(q, a)a(q) (163) 

c (2d + 2 7 (q2 + ^a2)) 

This equation is the nonlocal and nonlinear generalization of London equation in a long wavelength limit q — > 0, 
because the kernel Q(q, a) is function of q and magnetic field a(q). However the equation (|163p is correct when a 
magnetic field is much weaker than the critical field H <C He- 

Let's neglect by the nonlinearity, that is we suppose that the kernel Q is function of q only. Then we have: 

p 2 h 2 (T) p 2 A 2 (T) p 2 

j(q) = -27-— 4r^a(q) = -2 7 - - ^ \ 2 a(q) « -2 7 -A 2 (T) (l - / 2 q 2 ) a(q) = Q(?)a(q), (164) 



c (2d + 2 7 q 2 ) 



f q 2 

do 1 



where we took into account a slowness of changes of a gap in space: l$q -C 1. 2^- = 2 l °"jJ 2 — I 2 is a coherent 
length at temperature T = 0. Thus we /iawe t/ie nonlocal kernel Q(q), where radius of a nonlocality is equal to the 
coherent length Iq. This result corresponds to nonlocal Pippard electrodynamics (long wavelength limit). This fact 
proves nonlocality of the obtained functional of free energy of a superconductor (|141[ [T43|) . For generalization in case 
of a large value loq it is necessary to expand the free energy (|141|) in degrees of q. Then we can obtain a short 
wavelength limit of Q: Q ~ l/qlo- Starting from correctness of the asymptotics (|145[) and (|155j) of the functional 
(|141[ 1143ft we can make a conclusion about its correctness for description of a superconductive phase. 



IX. CONCLUSION. 



In this paper on the example of superconductivity we described the type II phase transition on a microscopic 
level, namely starting from first principles. This means, that the method of calculation of a free energy f2(T, N/V)) 
has been developed in a range of temperatures, which includes a point of pase transition, without introducing any 
artificial parameters of type of order parameter and sourscs of ordering, but starting from microscopic parameters of 
Hamiltonian only. Moreover, the theorems about connection of a vacuum amplitude with thermodynamics potentials 
are realized. 

Microscopic picture of a phase transition lies in the following. At switching of attraction between particles of a Fermy 
system the instability relatively formation of bound states of two fermions rises. The given states are characterized 
by Bethe-Solpiter amplitudes - amplitudes of pairing, which can be found exactly in the case of an isolated pair. 
However, in consequence of statistical correlations between pairs the amplitudes are determined by dynamics of all 
particles of a system. Their observing value is result of averaging over the system. Thus, a collective (condensate) 
of pairs exists. A particle propagating through a system interacts with fluctuations of pairing. As a result of such 
interaction a dispersion law of quasi-particles is changed and anomalous propagators appear. This means, that a 
spontaneous symmetry breakdown takes place. After consideration of interaction of particles with fluctuations of 
pairing all characteristics of a system must be calculated over the new vacuum with broken symmetry. So, calculation 
of a vacuum amplitude over new ground state gives a possibility to use the theorem about connection of a vacuum 
amplitude with a ground state energy (with free energy at nonzero temperature). As a result, the free energy is 
function of amplitudes D, = J7(AA + ), and their observed value minimizes the free energy. Analyze of the obtained 
functional of free energy shows, that the amplitude of pairing plays a part of an order parameter. Namely, the order 
parameter is Bethe-Solpiter amplitude averaged over a system due to statistical and dynamical correlations. Since in 
Nambu-Gor'kov formalism (the method of anomalous propagators) any phase transition can be described then 
our method can be generalized to the rest transitions (ferromagnetism and antiferromagnetism, waves of charge and 
spin density, ferroelectricity and so on). 



30 



The functional of a superconductor's free energy (|141[ 1143[) has been obtained in this paper using the developed 
method of microscopic description of phase transitions and generalizing its in the cases of spatial inhomogeneity and 
presence of magnetic field. The functional generalizes Ginzburg-Landau functional for cases of arbitrary temperatures, 
arbitrary spatial inhomogeneities and a nonlocality of a magnetic response. The equations of superconductor's state 
are extremals of the functional obtained by variation over the gap A and the magnetic field a. The equations 
determined equilibrium configurations of a gap and a magnetic field at given conditions. 

Appendix A: The method of uncoupling of correlations and Dyson equation. 

As it was shown in [20J for propagators G (one-particle), K 2 (two-particle), K 3 (tree-particle) ... K n (n-particle) 
the set of coupling equations can be written. The set of equations is analogous to BBGKY hierarchy for a s-particle 
probability density function. For example, for G and K 2 these equations have a view: 



G(l,l') = G (l,l')+W G (l,2)C/(2, 3)^(2, 3; l',3 + )dx 2 dx 3 (Al) 
A" 2 (l, 2; 1', 2') = G(2, 2')G (1, l') - G(2, l')G (l, 2') + * f G (l, 3)C/(3, 4)^(3, 2, 4; 1', 2', A+)dx 3 dx 4 , (A2) 



where 3+ = (£3, t 3 + 0), [7(2, 3) = V(&, ^)5(t 2 - t 3 ), dx = d£dt. The cross term G(2, 2')G (1, 1') - G(2, l')G (l, 2') 
appeared as result of calculation of Fermy symmetry of particles (it would be " + " for bosons). The equation for a 
three-particle propagator K 3 will be is determined by four-particle propagator K4 and so on: K n — f(K n +i)- 

It is obviously that this set of equations can not be solved. However, in most cases for description of a system it is 
enough to know functions G and K 2 (less). Then the method of uncoupling of correlations is used. The function K 2 
in a zero approximation is 

K$ ] = G (l, l')G (2, 2') - G (l, 2')G (2, l'). (A3) 

We can see, that in an absence of interaction the two-particle propagator is represented in a multiplicative form by the 
one-particle free propagators Go- Statistical correlation exists only in the course of Pauli exclusion principle. Then 
in first approximation let's use the free propagator K 2 °^ instead of the dressed propagators K 2 in the formula (|A1[) . 
Hence the correction for the dressed one-particle propagator G is 

G (1) (l,l') = % { Go(l,2)J7(2,3)Go(3,3 + )G (2,l')dz; 2 cte3-i f G (l, 2)17(2, 3)G (3, l')G (2, 3)dx 2 dx 3 . (A4) 



The correction (|A4I) is represented graphically in FigfTOJ Integration is carried over coordinates of the internal lines. 
The first term corresponds to direct Hartree interaction, the second term corresponds to exchange Fock interaction. 
In order to obtain the next approximation for G it is necessary to find In a symbolic representation we have 
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Figure 10: The correction of first order Gi obtained by uncoupling of correlations in the equation (|A1[) . 
the equations: 

K { 2 ] = G (1) G - G G (1) + iG UK { 3 0) 

G (2) =iG UK { 2 \ (A5) 

The procedure of uncoupling of correlations can be represented in another way. Let a correction for two-particle 
propagator is determined by the matrix elements: Vkiu is a direct interaction and Vuik is an exchange interaction. 
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A two-particle propagator has two entering momentums (represented by lines with corresponding indexes) and two 
outgoing momentums. An one-particle propagator has one entering momentum and one outgoing momentum. The 
procedure of uncoupling of correlations consist in the fact, that we connect two lines in K taking into account 
conservation of momentum and spin: Fig llll The connection means integration over intermediate momentums and 
energy parameters (in a momentum-energy representation £,i — > k, ui) in the formula (|A4[) . As a result, we have the 
same diagrams for G as in FigfTOl 






Figure 11: The graphical method of uncoupling of correlations. In the diagram for the first order correction K\ of two-particle 
propagator we connect two lines taking into account conservation of momentum and spin. As a result, we have the first order 
correction Gi for one-particle propagator. 

The procedure can be generalized to higher corrections (with two and more lines of interaction). In general case 
the rules of diagram technique are: 

1. The multiplier iG(k, lo) is associated with each bold line, the multiplier zGo(k, cj) is associated with each thin 
line 

2. The multiplier —iVkimn is associated with each dashed line of interaction between particles, and the multiplier 
—iVki is associated with line of interaction of a particle with an external field. 

3. The multiplier —1 associated with each fermion loop and we make summation over all possible spin configura- 
tions. 

4. Momentum k + 1 = m + n, energy parameter and spin are reserved in every vertex. 

5. Integration is made over each intermediate momentum and summation is made over each intermediate energy 
parameter: — > y J {V is volume of a system) and J 

The diagrams of the type Fig llOl are summarized with help of the mass operator £ - any diagram without external 
lines. Hence a dressed propagator G can be found from Dyson equation: 

iG = iG + iG (-i)EiG ^> G= — (A6) 

G - £ 

The mass operator £ has the sense of a mean field of all particles of a system acting on a test (marked) article. In 
this fact the sense of the procedure of uncoupling of correlations is: interaction and propagation of all particle of a 
system is reduced to propagation of each particle in the mean field of all rest particles. 

Other approach exists (more widely represented in literature) for obtaining of the diagram expansion for G. In this 
approach a one-particle propagator is determined as 

-i{^\T[U{T 2 ,T 1 )C^{t 2 )C+ (t 2 )]|$ ) 

G(k 2 ,ki,t 2 -ti)= lim = £2 , (A7) 

T x -oo(l - iS) <*o|J7(22,Ti)|*o> 
T 2 -> +oo(l - id) 

where C + (t),C(t) is creation and annihilation operators in interaction representation, U is evolution operator in 
interaction representation, $ is ground state of a system of noninteracting fermions. This definition of a propagator 
is equivalent to the definition 

G(k,t a -t 1 ) = -i<* |T[O ki(r (t a )Cj |(r (t 1 )]|*o>, (A8) 

where C + (t),C(t) is creation and annihilation operators in Heisenberg representation, <F is ground state of a system 
of interacting fermions, and expansion of G in series of perturbation theory is possible if the condition of adiabaticity 
is realized (|8]): 

($o|*o}^0. (A9) 
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In the method of uncoupling of correlations we didn't use the condition (|A9|) and Wick theorem unlike the standard 
formulation of the perturbation theory based on (|A7[) . Thus the advantage of stated above method of uncoupling of 
correlations consists in that we can formulate a perturbation theory without using of the adiabatic hypothesis (|A9j) . 



Appendix B: The method of uncoupling of correlations for a vacuum amplitude in the case of normal 

processes. 

Let us consider the processes of direct and exchange interaction of first order, which is described by the matrix 
elements Vkm and Viku > moreover the interaction doesn't act to spins of particles. Contribution to vacuum amplitude 
of such processes is (let's suppose t-x > t\ for dchniteness) : 

R(t) = 1 + ^7 / (-^v klkl ) (^oiq*- (t 1 )c+ cr (t 1 )q kjQ1 (t 1 )q ljjS (t 1 )|*o) 

' a >P k,l 

+Tf f **EE (-l v ^i) (*o|cJ a (ii)cf a (t 1 )c k)a (* 1 )q ia (t 1 )|$ ) 
' J0 « k,l v 7 

+ ijf* J fAEE(-5««)EE(-5'4'«" 

,y0 J0 a,/3 k,I V 7 «'./8'k',l' 

4/^/^ee(-^)ee(-^ 



Q kj x 7 «' k',1' 



2! 



2 / ^EE(-^-)EE (-^w ) 



a,/3 k,l V 7 <*' k',1' 

x <*o|C+, (taJCf a X*a)^^(t 2 )C r ^(t a )Cf fl (t 1 )C+ a (t 1 )C kia (t 1 )q i/J (i 1 )|*o) 



(Bl) 



For approximate calculation R(t) we shall use the method of uncoupling of correlations, which lies in the fact that an 
average of four creation and annihilation operators is represented by a product of averages of pairs of the operators: 
(C + C + CC) — > (C + C)(C + C). The averages correspond to propagators of particles with initial and final states 
corresponding to the matrix element of interaction Vkimn taking into account conservation of momentum and spin. It is 
achieved by preliminary transposition of the operators before the uncoupling taking into account Fermy commutation. 
For Hartree and Fock processes the procedure corresponds to the diagrams in Fig |12l In each process of scattering we 
connect incoming and outgoing lines taking into account the laws of conservation. The obtained diagram must not 
have free ends. Then we have the amplitude of transition " vacuum- vacuum" . Analytically it will be so: 



i?(t)«i + (-i) 2 i f^EE(i^) (*o|c , + /3 (i 1 )q /; (t 1 )|<i> )($o|C'+ Q (t 1 )c k!Q (i 1 )|<i>o) 

• J ° a ,f3 k 1 V / 

/VEE (-> ifeW ) ^l c, L( <1 ) c k, Q (*i)l $ o>(<i>o|c 1 + a (< 1 )q a (i 1 )|<i>o) 
' Jo " k,l V J 

+(-d 4 ^ fdt 2 T^EE (-M E E (-W 

Z - J ° J ° a ,P k,l V J «*',/>' k',1' V 1 J 

A*Q\<$fil)Cl tP (t x )\* )(* \^ 

+{ - lf h.[ * jf E E (->) E E (-| w ) 

JO Jo a k l \ / a , k , j, \ / 
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2!< 



7 , \ ^ \ ^ / t T r \ \ ^ \ ^ / t T r 

/3 k,i v " k r 1 

+ . . . = 1 + (i?f artree + R{ ock ) + (i?f artree + i?f ocfc ) 2 + . . . = exp(i2f ar * ree + i?f ocfc ) (B2) 
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Figure 12: The procedure of uncoupling of correlations in vacuum amplitude for a correction of first order. In the average of 
matrix element of interaction operator V(C + C + CC) we connect the lines taking into account conservation of momentum and 
spin. As a result, we obtain the correction of first order Ri for a vacuum amplitude of a view V (C + C) (C + C) . 

The uncoupling lets to write expression for R\{t) via normal propagators: 



Hartree 



(2s + 1 



) 2 /**e 
j ° k.i 



--jVklkl 



iG a (\,t - t)iG (\s.,t - t) 



{-2i)^V km B {l)B {k)t 
k.l 



(B3) 



R{ 



Fock 



(2* + l)(-l) / dtJ2 

J ° k,l 

iJ2v im Bo(l)B (k)t. 
k,l 



-7,Vi kld ) iG (l,t-t)iG (k,t-t) 



(B4) 



The multiplier (2s + 1) is result of summation over spin states (number of spin configurations), s — 1/2. We can see, 
that the method of uncoupling of correlations lets to calculate vacuum amplitude simply, selecting contributions of 
processes of each type. The method can be generalized to the processes of higher order. We can see, that the proposed 
method of obtaining of diagram expansion for a vacuum amplitude by uncoupling of correlations doesn't demand of 
an use of Wick theorem and the adiabatic hypothesis. 
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